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Abstract 

We investigate the influence of hadron injection from evaporating primor- 
dial black holes (PBHs) in the early stage of the primordial nucleosynthesis 
era (t ~ 10^'^ — 10^ sec). The emitted quark-antiquark pairs or gluons im- 
mediately fragment into a lot of hadrons and scatter off the thermal plasma 
which is constituted by photons, electrons and nucleons. For the relatively 
low mass holes we point out that the dominant effect is the inter-conversion 
between ambient proton and neutron through the strong interaction induced 
by the emitted hadrons. Even after the freeze-out time of the week inter- 
actions between neutron and proton, more neutrons are produced and the 
synthesized light element abundances could be drastically changed. Com- 
paring the theoretical predictions with the observational data, we constrain 
the PBH's density and their lifetime. We obtain the upper bound for PBH's 
initial mass fraction, /3 < lO^^o f^j. iqS^ < M < lO^^g, and /? < 10~22 for 
lO^Og < M < 3 X lO^Og. 
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I. INTRODUCTION 



Primordial black holes (PBHs) are formed in the hot early Universe if overdensity of order 
of unity exists and a perturbed region enters the Hubble radius . They serve as a unique 
probe of primordial density fluctuations on small scales. For example, if massive compact 
halo objects (MACHOs) found by gravitational microlensing observation toward the 
Large Mageranic Clouds turn out to be the PBHs, which originate from density fluctuations 
generated during inflation, we can determine model parameters of the inflation model with 
high accuracy 0. In fact, this possibility is attracting more attention these days because 
none of the conventional candidates of MACHOs are plausible from various astrophysical 
considerations and observations. 

Even the nonexistence of PBHs over some mass ranges, however, would provide useful 
informations on the primordial spectrum of density fluctuations. In this sense it is very 
important to obtain accurate constraints on the abundance of PBHs on each mass scale. 

Cosmological constraints on the mass spectrum of the PBHs can be classified to three 
classes. The first one applies to heavy black holes with mass M > 4 x 10^'^g which have 
not evaporated by now. Current mass density of such holes should not exceed the total 
mass density of the universe. More stringent constraints are obtained from the microlensing 
experiments on holes with sub-solar masses 0. The second class is due to the radiation 
of high energy particles from evaporating black holes Various constraints have been 

imposed from primordial big-bang nucleosynthesis (BBN) |§-[TD|, microwave background 
radiation [jll| , and gamma-ray background radiation . Finally, there may be yet another 
class of constraints if PBHs do not evaporate completely but leave relics with mass of order 
of the Planck mass or larger |T^. Their mass density should remain small enough. 



The purpose of this paper is to reanalyze the effects of evaporating primordial black holes 
on BBN in order to improve constraints on their mass spectrum. Although extensive work 
was done on this issue in 1970s just after the idea of PBH was proposed, to our knowledge, 
it has not been studied almost two decades now, and old constraints are still used in the 
literature |14,|15[. In order to compare the previous work with the modern view of BBN and 
high energy physics, we start with a brief review of what has been done on the subject. 

Vainer and Nasel'skii studied the effects of the injection of high-energy neutrinos 
and antineutrinos which change the epoch of the freeze out of the weak interactions and 
also the neutron-to-proton ratio at the onset of the nucleosynthesis. This results in the 
increase of ^He. Demanding that its primordial abundance, Yp, should satisfy Yp < 0.33, 
they concluded the ratio of the energy density of PBHs to that of baryons should be smaller 
than 10 — 10^ for holes with mass M = 10^ — 3 x lO^^g. This approximately corresponds to 
/3{M) < 10^22 - 10"^^, where (3{M) is the initial fraction of PBHs with mass M to the total 
energy density of the universe when the horizon mass is equal to M, i.e. [3{M) = {pbu/ Ptot)i- 
These authors also studied the effects of entropy generation from PBHs but they obtained a 
modest constraint that the energy density of the relevant holes should be smaller than that 
of the photons at evaporation. 

More detailed numerical analysis of the latter effect was done by Miyama and Sato [0, 
who calculated the effect of entropy production from PBHs with mass M = 10^ — lO^^g 
evaporating during or after the nucleosynthesis. If these PBHs were too abundant, the 
baryon-to-entropy ratio at the nucleosynthesis should be increased which would result in 
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overproduction of ^He and underproduction of D. They demanded that Yp < 0.29 and that 
mass fraction of D should be larger than 1 x 10~^. If we approximate their result with a 
single power law, we find 

f3{M) < 10-^^M^o^/^ (1) 

for M = 10^ - lO^^g, where Mio = M/IO^^ g. 

Zel'dovich et al |^ , on the other hand, studied the effect of emission of high-energy nu- 
cleons and antinucleons from PBHs. They claimed such emission would increase deuterium 
abundance due to capture of free neutrons by protons and spallation of "^He by emitted 
particles. Their conclusions were 

10^°g, 
5 X 10^°g, 

1° - 5 X 10"g, (2) 
lO^^g. 

Vainer, Dryzhakova, and Nasel'skii then performed numerical integration of nucle- 
osynthesis network in the presence of neutron-antineutron injection, taking the neutron 
lifetime to be 918(±14) sec and using the observational data Yp = 0.29 ± 0.04 and the mass 
fraction of D to be 5 x 10~^. They calculated spallation of ^He and resultant extra produc- 
tion of D assuming that PBHs were produced out of scale-invariant density fluctuations. As 
a result they obtained the constraint (3 < 10~^^ for the relevant mass range. 

Lindley [|1^] considered another effect of evaporating PBHs with mass M > lO^^g, namely, 
photodissociation of deuterons produced in the nucleosynthesis. He found the constraint 
/? ^ 3 X 10~^°Mio^, and concluded that photodestruction was comparable to the extra 
production of deuterons discussed in §]. 

Since these papers were written, both observational data of light elements and the neu- 
tron lifetime have changed considerably [|T6|. More important, however, is the change in our 
view of high energy physics of the energy scale relevant to PBHs evaporating in the nucle- 
osynthesis era. That is, hadrons are not emitted in the form of nucleons or mesons but as a 
quark-gluon jet in the modern view of quantum chromodynamics (QCD). Thus we should 
adopt the Elementary Particle Picture of Carr ||^ and assume that elementary particles in 
the standard model are emitted from PBHs and they generate jets. In fact, the number of 
jet-generated particles far exceeds that of directly emitted counterparts. In this sense the 
previous calculation of BBN in the presence of evaporating PBHs should be entirely revised. 
The above procedure has already been taken in the analysis of cosmic rays from evaporating 
PBHs |1^. We incorporate emission of various hadrons from PBH-originated jets to BBN 



for the first time. 

We adopt a simple and conventional view that PBHs are produced with a single mass, M, 
when the horizon mass is equal to M, namely at t = 2.4 x 10~^^Miosec = tform, and obtain 
an improved constraints on their initial abundance at each mass scale. Recent numerical 
calculations of PBH formation, however, have revealed the mass spectrum of PBHs spreads 
rather widely even if the spectrum of primordial density fluctuation is sharply peaked on a 
specific scale |1^JT^. In particular, the authors of [|18| discovered the critical behavior and 



its cosmological consequences have been discussed in [^. By convolving our results with 
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the mass functions obtained in these papers, one can obtain constraints on PBH formation 
in more reahstic situations. 

When PBHs evaporate and emit various particles in BBN epoch, i.e. at T = lOMeV — 
IkeV, such high energy particles interact with the ambient photons, electrons and nucleons 
and they finally transfer all the kinetic energy into the thermal bath through the electro- 
magnetic interaction or the strong interaction. Through the above process, such dangerous 
high energy particles possibly induce the various effects on the background and change the 
standard scenario considerably. 

Once quark- ant iquark pairs or gluons are emitted from a PBH, a lot of hadrons, e.g. 
pions, kaons and nucleons (protons and neutrons) are produced through their hadronic 
fragmentation. They inter-convert the ambient protons and neutrons each other through 
the strong interaction. If the inter-conversion rate between neutrons and protons becomes 
large again after the freeze-out time of the weak interactions in the standard BBN (SBBN) 
or t ~ 1 sec, the protons which are more abundant than the neutrons at that time are 
changed into neutrons. That is, this results in an excess of neutrons compared with SBBN. 
Therefore, the hadron injection significantly influences the freeze-out value of the neutron- 
to-proton ratio and the final abundances of ^He, D and ^Li are drastically changed. In this 
case it is expected that both ^He and D tend to become more abundant than in SBBN. If 
PBHs are so massive that they continue to emit particles even after "^He are produced, then 
spallation of ^He due to high energy particles will tend to increase the final D [} In this 
paper, as a first step of the full analysis, we concentrate on the low-mass PBHs and consider 
only the former effects. The latter issue will be addressed in a separate publication . 



Reno and Seckel |^ investigated the detail of the physical mechanism and the influences 
of the hadron injection from long-lived massive decaying particles on BBN. They constrained 
the parent particle's lifetime and the number density comparing the theoretical prediction 
of the light element abundances with the observational data. Here we basically follow their 
treatment and apply it to the hadron injection originated in the PBH evaporation. 

The rest of the paper is organized as follows. In §11, we review properties of black hole 
evaporation and jets. In §111 formulation to calculate the effect of hadron injection to BBN 
is given and observational data of light element abundances are summarized. The result is 
presented in §IV. Finally §V is devoted to the conclusion. We use the units c = h = ks = 1. 

II. EVAPORATION AND JETS 

First we briefly summarize basic results of PBH evaporation. As was first shown by 
Hawking a black hole with mass M emits thermal radiation with the temperature given 
by 

More precisely, a neutral non-rotating black hole emits particles with energy between Q and 
Q + dQ at a rate 



^This situation has been studied in |24| in the context of late-decaying particles. 
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where s is the spin of the emitted particle. Tg is its dimensionless absorption coefficient 
whose functional shape is found in |^T|. It is related with the absorption cross section 
as{M,Q) as Ts{M,Q) = Q'^as{M,Q) /it. In the high-energy limit Q^T, ag approaches to 
the geometric optics limit ag = 277rG^M^. 

The average energies of a neutrino, an electron, and a photon are given by E,y = 4.22Tbh, 
Ee = 4.18Tbh, and = 5.71Tbh, respectively. The peak energy of the flux and that of the 
power are within 7% of the above values |^ . Averaging over degrees of freedoms of quarks 



and gluons in the standard model, we find the average energy E = 4.4Tbh- 

MacGibbon elaborated on the lifetime of a PBH, tbh, summing up the contribution of 
all the emitted particles and integrating the mass loss rate 

= -5.34 X 10-V(M)Mro' sec-\ (5) 



over the lifetime ||2^. The result is approximately given by 

tbh = 6.24 X 10V(M)-^Mj^0 sec. (6) 

Here /(M) is a function of the number of emitted particle species normalized to unity for 
M ^ lO^^g holes which emit only photons, three generations of neutrinos and anti-neutrinos, 
and gravitons. Depending on the spin s, each relativistic degree of freedom contributes to 
/(M) as g 



/,=o = 0.267, = 0.060, /,=3/2 = 0.020, /,=2 = 0.007, (7) 

/^^-^/2 = 0.147 (neutral), /^=i/2 = 0.142 (charge = ±e). 

Summing up contributions of all the particles in the standard model, we find /(M) = 14.34. 
In this case the lifetime reads 

TBH = 435Mj^oSec, (8) 

or 

M = 1.32x lO^f^Vg. (9) 

VI sec/ 

Thus a PBH evaporating by the end of BBN t ~ 10^ sec has a mass M ^ 10^°g and the 
temperature Tbh ~ ITeV. 

For such a PBH with temperature higher than the QCD scale, Aqcd ~ 10^'^MeV, it has 



been argued that particles radiated from it can be regarded as asymptotically free at 
emission. The emitted quarks and gluons fragment into further quarks and gluons until they 
cluster into the observable hadrons when they have traveled a distance Aq^j-, ~ 10~^^cm. 
The hadron jet thus produced would be similar to that produced in e~^e~ annihilation. 
We now estimate the average number of the emitted hadron species Hi per jet as 



iV^' = /H.^, (10) 
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where (Nch) is the averaged charged-particle multiphcity which represents the total number 
of the charged particles emitted per two hadron jets, fn, is the number fraction of the hadron 
species Hi to all the emitted charged particles. 

It is reasonable to assume that the averaged charged-particle multiplicity is independent 
of the source of the hadron jet. In this paper we use the data obtained by the e^e~ collider 
experiments and extrapolate it to higher energy scales. A number of experimental data 
are available at least up to ^/s ~ 100 GeV where ^/s denotes the center of mass energy. 
Recently LEPII experiments (ALEPH, DELPHI, L3 and OPAL) yielded useful data for y/s 



= 130 — 172 GeV. One can fit the data by a function, {Nch) = a + bexp{c^^\n{s/ A'^)), where 
A is the cut-off parameter in the perturbative calculations associated with the onset of the 
hadronization and a, b and c are constants ^7^. The above functional shape is motivated 
by the next-to-leading order perturbative QCD calculations. We fit the data of the e~^e~ 
collider experiments for ^/s = 1.4 — 172 GeV P3] and we find 



(N^h) = 1-73 + 0.268 exp (lA2^\n{s/A^)) , (11) 



with = 90.3 for 93 data points and the error of the fitting is about 10%. Here we have 
taken A = 1 GeV. We assigned systematic errors of 10% for y/s = 1.4 — 7.8 GeV because 



MARK I and 772 experiments do not include the systematic errors in their data, see p8 
In Fig. [| we plot (Nch) for the center of mass energy ^/s. When we apply the results of the 
e~^e~ collider experiments to the case of PBH evaporation, we take ^/s = 2E as the energy 
of two hadron jets from the PBH evaporation because we assume that two hadron jets are 
induced per one quark-antiquark pair emission here. Then we can approximately estimate 
the average energy of an emitted hadron species Hi as Eh^ — E/ {Nch)- 

The next issue is to find fui- In this paper we adopt the experimental data of at 
y/s = 91.2 GeV, which is the highest energy for which these data are available |]28|, and we 



assume that do not change significantly in the energy range ^/s ~ 100 GeV - 20 TeV. [| 
From the table of //f. in we must estimate effective /^^ for each hadron species on the 
relevant time scale. That is, we must also take into account the decay products of those 
particles whose lifetime is too short to affect BBN. 

Let us estimate the time scale at which we should estimate The emitted hadrons 
do not scatter off the ambient nucleons directly. At first the emitted high energy hadrons 
scatter off the background photons and electrons because they are much more abundant 
than the nucleons. For the most part of the epoch at the cosmic time t ^ 10^ sec, as we can 
see later, it is expected that the emitted particles are quickly thermalized and they have the 
kinetic equilibrium distributions before they interact with the ambient nucleons. In addition 
as we show in Sec. it is reasonable to treat the emitted hadrons to be homogeneously 
distributed. Then we use the thermally averaged cross sections {crv)^'__^j^, for the strong 



^The above energy range is approximately the PBH's temperature corresponding to the lifetime 
tbh = 10~^ - 10^ sec. For ^/s > 350 GeV, t-t pairs would be produced and they would change 
the form of the charged particle multiplicity and the hadron fraction. Since we do not have the 
experimental data for such high energy regions, we extrapolate (N^h) to the higher energy regions 
and we take //f. as a constant. 
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interaction between hadron Hi and the ambient nucleon N, where N denotes proton p or 
neutron n. For a hadron interaction process N + Hi ^ N' + ■ ■ ■, the strong interaction rate 
is estimated by 

where n^v is the number density of the nucleon species A^, r]i is the initial baryon to photon 
ratio {= ns/n^ a.t T ^ lOMeV), ub = rip + Un denotes the baryon number density, /at = 
nN/^B, and Ti, is the neutrino temperature 0. 

Thus all we need to consider are particles with lifetime larger than O{10~^) sec. The 
corresponding mesons are vr"*", vr", K^, K~ , and Kl and the baryons are p, p, n, and n. 
We have therefore calculated the final yield of these particles per two jets, n^', out of Table 



38.1 of |28[. The results are: 

n^"" = 14.1, n^" = 14.1, 
n^^ = 1.67, = 1.67, n^^ = 1.19 (13) 

nP = nP = 0.772, n" = n" = 0.772. 

In estimating n" we have assumed that the number of neutrons emitted directly from a jet 
is the same as that of protons. Then we obtain the number fraction //f. by 

^""^ ^ (iV,,(v/i = 91.2GeV))- ^^^^ 



III. HADRON INJECTION AND BBN 

A. Hadron scattering off the background particles 

When there are sufficient electrons and positrons in the universe (T ^ 0.025 MeV), it is 
expected that such charged mesons are quickly thermalized through the Coulomb scattering. 
Since the stopping time to lose the relativistic energy is estimated as 

T^h ^ 10~^^ sec{E/ GeV)/(T/ MeV)^ (15) 

where E is the kinetic energy of a charged meson, the long-lived charged mesons are ther- 
malized and scatter off the ambient nucleons by the threshold cross section before they decay 
The thermally averaged cross sections for vr^ are obtained by [EH 



^The Ti, dependence comes from that the baryon number density decrease as oc R{t)^^ and 
the scale factor exactly decrease as R{t) oc while the universe adiabatically expands. If large 
entropy is produced, the above relation must be changed. In the hadron injection scenario, however, 
the large entropy production region is severely excluded. Therefore Eq. ([l^ ) is a good estimate of 
the time scale of the hadron interactions. 
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{av)ll^ = 1.7 mb, (16) 
{av);-^^ = 1.5C-(T) mb, (17) 

where C^'{T) is the Coulomb correction factor. Because the reaction p+ + tt^ n + ■ ■ ■ is 
enhanced due to the opposite-sign charge of the initial state particles, we correct the strong 
interaction rates by multiplying C^* (T) to that which are obtained by ignoring the Coulomb 
corrections. The Coulomb correction factor is estimated by 



where ^j(T) = ay fii/2T, a is the fine structure constant and fii is the reduced mass of the 
hadron Hi and the nucleon. 

The cross sections for K~ are obtained by p6| 

(19) 

(20) 

(21) 
(22) 

Following Reno and Seckel ^6] we ignore interaction because n + K'^ p + is the 





= 26 mb. 




= 34 mb. 




= 31C^(T) mb. 


Wv)p^p 


= 14.5C^(T) mb 



endothermic reaction which has Q = 2.8 MeV. 

On the other hand, among neutral kaons, Kl has a long lifetime O{10~^) sec. Since 
Kl does not stop through the electro-magnetic interactions, we adopt the following strong 
interaction cross sections which are obtained by the initial energy distribution in the hadron 



fragmentation [^] 



(^^)nip = 7 mb, (23) 

{av)i^^^ = 10 mb, (24) 

{avYX = 7 mb, (25) 

{(xv)^^p = 10 mb. (26) 



As for the emitted high energy nucleons, we should treat the thermalization process more 
carefully because the stopping process is not so simple. Proton and antiproton are stable and 
we should worry about the efficiency to lose the kinetic energy at the lower temperature. 
At least T ^ 0.02 MeV, protons are quickly thermalized through Coulomb scattering off 
the electrons and Inverse Compton scattering off the photons. For t ^ 3 x 10^ sec, on the 
other hand, the stopping process of protons proceeds through the nuclear collisions with the 
ambient protons and light elements. In this case such high energy protons may induce the 
"^He fissions because ^He has already been synthesized at around t ~ 300 sec. For neutron 
and antineutron the efficiency of the thermalization is more severe. Neutrons are efficiently 
stopped by the electron scattering until T ~ 0.09MeV. Thus the late time emission of the 
high energy neutron would induce the light element fissions. However, in this paper we treat 
the neutrons as if they are approximately thermalized in the entire parameter range because 
we are primarily concerned with the effects of low-mass PBHs here. We follow the Reno 
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and Seckel's treatment that a nucleon-antinucleon pair is regarded as a meson NN. Then 
the NN meson induces the inter-converting reactions hke + NN N' + ■ ■ ■ and the 
thermally averaged cross sections are given by Ef 



{av):tn = 37 mb, (27) 

{av);t^ = 28 mb, (28) 

{av)^.^^ = 28 mb, (29) 

{(^v)ltp = 37 mb. (30) 

The above treatment may underestimate the deuterium abundance because it will be pro- 
duced by the hadro-dissociation of ^He if PBHs are so massive that they continue to emit 
high-energy hadrons even after ^He are formed. The hadron-induced dissociation process of 



the light elements will be discussed in the next paper |2^. For the other mesons and baryons 
whose lifetimes are much shorter, e.g. tt" —>■ 27 with r^o ~ O{10^^^) sec, they quickly decay 
into the standard particles and do not influence the standard scenario. 

On the other hand, for the even longer lifetime tbh ~ 10^ sec, there is another in- 
teresting effects on BBN. The emitted photons or charged particles induce the electro- 
magnetic cascade showers and produce many soft photons. Their spectrum has a cutoff at 
^max ^ mg/(22T), where mg is the electron mass ^9]. If E^^^ exceeds the binding energies 



of the light elements, these photons dissociate light elements and change their abundances. 
In fact, the energy of the photon spectrum which are produced by the PBH evaporation at 
t ^ lO'^(lO^) sec exceeds the deuterium ("^He) binding energy B2 — 2.2 {B^ ~ 20) MeV. In 
this case PBH's number density and the lifetime are severely constrained by the observa- 
tional data of the light element abundances JTU . 



B. Formulation 



As we mentioned in the previous section, the hadron emission mechanism by the PBH 
evaporation at t ^ 10^ sec induces extra interactions between emitted hadrons and ambient 
nucleons. That is, they enhance the inter-converting interaction rates between neutron and 
proton even after the weak interactions has already frozen out in the standard scenario and 
the freeze out values oi n/p ratio can be increased. The time evolution of the PBH's mass 
is given by 



M{t) 







(for 
(for 



t < 



tbhJ 
< 



Tbh ^ t) 



where M is the initial mass of PBH when it was formed. 

Then the time evolution equations for the number density of a nucleon A^(= 
represented by 



dn 



N 



dt 



+ 3H{t)nN 



dn 



N 



dt 



K 



(31) 



p, n) is 



(32) 



J SBBN 



where H{t) is the cosmic expansion rate, [driN / dt]sBBN denotes the contribution from the 



standard weak interaction rates []I6[ and nuclear interaction rates, is the hadronic branch- 
ing ratio, J{t) denotes the emission rate of the hadron jet per unit time and K^-^n' denotes 
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the average number of the transition N ^ N' per one hadron jet emission. The emission 
rate of the hadron jet is estimated by 

^ ' E{t) dt ' ^ ' 

where nBH(^) is the number density of the PBHs. 

Though PBHs generally emit not only quarks and gluons but also all lighter particle 
species than the temperature of PBH, the emitted neutrinos, photons and the other electro- 
magnetic particles do not influence the light element abundances significantly for the rela- 
tively short lifetime. As we noted in the previous section, for tbh ~ 10^ sec the injection 
of photon or the other electro-magnetic particles do not induce the photo-dissociation of 
the light elements. On the other hand, the emitted neutrinos scatter off the background 
neutrinos and produce the electron-positron pairs. Although they also induce the electro- 
magnetic cascade, we should not be worried about photo-dissociation for the same reason. 
Hence we concentrate on the effects of hadron injection in BBN epoch. Then the resultant 
upper bound on the abundance of PBHs, (3, turns out to be proportional to Bf^^ . Below we 
analyze the case Bh = 1, because its magnitude is not precisely known, although we expect 
Bh = 0(0.1 — 1). Anyway the constraint for other cases with Bh ^ 1 can easily be obtained 
from the above scaling law. 

The average number of the transition N ^ N' per jet is expressed by 

K^^N' = J2n'''Rn^n', (34) 

where runs the hadron species which are relevant to the nucleon inter-converting reac- 
tions, N^' denotes the average number of the emitted hadron species Hi per jet which is 
given by Eq. (|l^) and R^^j^, denotes the probability that a hadron species Hi induces the 
nucleon transition N'. The transition probability is estimated by 

dec ~^ abs 



where P^^ = t^^ is the decay rate of the hadron Hi, P^^ is the total absorption rate of 
Hi. For Kl, which is not stopped, the decay rate is approximately estimated by P^^ = 
ttikl/ Eki^Tj^I where Ek^ is the averaged energy of the emitted Kl. 



C. Observational constraints 

As we mentioned in the previous subsection, the hadron injection increases the freeze-out 
value oi n/p ratio. The remaining neutrons are included in the deuteriums rapidly and it 
burns into ^He. Since the effects of the late time hadron emission tend to increase D and ^He 
we can constrain mass and abundance of PBHs comparing the yield with the observational 
light element abundances. In this subsection we briefly review the observational data of the 
light elements D, ^He and ^Li. 

The primordial deuterium is measured in the high redshift QSO absorption systems. 
Because it is expected that such a Lyman limit system is not contaminated by any galactic 



10 



or stellar chemical evolution, we regard it as a primordial component. At present we have 
two class of the observational D values, Low D and High D. However, since it is premature 
to determine which component is primordial, we consider both of two cases in this paper. 
Buries and Tytler observed clouds at z = 3.572 towards Q1937-1009 and a.t z = 2.504 towards 



Q 1009+2956 and they obtained pO] 



D/H = (3.39 ±0.25) X 10"^ : Low D, (36) 

where D/H denotes the deuterium number fraction and the error is the la value. On the 
other hand, Webb et al. reported the higher abundance in relatively low redshift absorption 
systems at ^ = 0.701 towards Q1718+4807, D/H ~ 0(10"^) |3l|. Tytler et al. also observed 
the same clouds and obtain the high deuterium abundance independently p2| . 



D/H = (0.8 - 5.7) X 10-^ : High D, (37) 

where the uncertainty is 2a. 

The primordial component of the '^He mass fraction Yp is measured in the low metalicity 
extragalactic HIT region. In addition to the primordial component, ^He is also produced in 
stars together with Oxygen and Nitrogen through the stellar evolution mechanism. Therefore 
in order to obtain the primordial component from the observational data we should regress 
to the zero metalicity for the observed ^He values. Olive, Skillman and Steigman ||3^ 



adopted 62 blue compact galaxies (BCG) observations and they obtained the relatively 
"Low" abundance, Yp ~ 0.234. However, Thuan and Izotov |Q pointed out that the effect 
of the Hel stellar absorption, which was not considered properly in |Q , is very important 



and they reported the relatively "High" value, Yp = 0.245 ± 0.004. Recently Fields and 
Olive reanalyzed the observational data and they reported 

Yp = 0.238 ± {0m2),tat ± (0.005),,,^, (38) 

where the errors are the la values. 

The primordial value of ^Li/H is measured in the Pop II old halo stars. Since the low 
mass stars , i.e. the low temperature stars, have the deep convective zone, the primordial 
component should be considerably depleted in the warm interiors. For the large mass stars 
which have the relatively high effective surface temperature Te// ^ 5500K, it is known that 
the primordial components are not changed and they have a "plateau" of the ''Li/H as a 
function of T^jf. On the other hand it is also reported that ^Li/H abundances decrease 
along with Iron fraction [Fe/H]. For the plateau stars, however, the ''Li abundances are 
not changed with decreasing [Fe/H] for [Fe/H]^ —1.5. Bonifacio and Molaro observed 41 
plateau stars and reported p6| 



logio('Li/H) = -9.76 ± {0.012)stat ± (0.05),,,^ ± (0.3),^^, (39) 

where ^Li/H denotes ^Li number fraction, the errors are 1 a values and we take an additional 
systematic error A logj^QC^Li/H)^^^ = 0.3, for fear that we may underestimate the stellar 
depletion and the production by the cosmic ray spallation. 

In order to get the conservative bound for the hadron injection induced by PBH evapo- 
ration, we adopt the following mild observational bounds here. 
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Yp < 0.252 (2(t), 

D/H < 4.0 X 10"^ (2ct) for LowD 
D/H < 5.7 X 10"^ (2(t) for HighD, 
3.3 X 10"" < ^Li/H < 9.2 x 10"^° (2a) 



(40) 
(41) 
(42) 
(43) 



where we summed all the errors in quadrature. 



IV. RESULTS 



In this section we compare the theoretical predictions of the light-element abundances 
in the hadron injection scenario with the observational constraints. Now we have three free 
parameters, the baryon to photon ratio r], the PBH's lifetime tbh, and the initial number 
density of the PBH normalized by the entropy density, s, Ybh = n-sn/s. We start the BBN 
calculation at the cosmic temperature T = 100 MeV. Since t] is the value at present time, 
the initial value rji should be set to an appropriate value which turns out to the present rj 
after the possible entropy production due to PBH evaporation and the photon heating due 
to e~^e~ annihilation. As we noted in the previous sections, the lifetime tbh characterizes not 
only the decay epoch, but also the PBH's initial mass M and the initial PBH temperature 
Tbh- Here we can relate Ybh to the initial mass fraction (3 by the following equation: 



Since the hadron injection tends to increase the produced D and "^He abundances in the 
present situation, the parameter range of t] is necessarily restricted to a narrow region if we 
adopt a set of the observational upper bounds for D and "^He. First we discuss the case of 
"LowD" and take Eqs. (|40|), (^l]) and (^3]) as a first set of the observational constraints. 
In this case the baryon-to-photon ratio is restricted in r] = (4.7 — 8.6) x 10~^°. 

In Fig. we plot the upper bounds for /5 which come from the observational constraints 
of ^He, D and ^Li for the PBH mass M = 10^ - 3 x lO^^g at t] = 6.0 x 10"^°. The above 
mass range corresponds to the lifetime tbh = 4.4 x 10~^ — 10^ sec. As we noted earlier, the 
shorter lifetime, tbh ~ 10"^ sec, do not affect the freeze-out value of n/p and do not change 
any predictions of SBBN. However, if the lifetime becomes tbh ^ 10~^ sec, the freeze-out 
value oi n/p ratio is increased by the hadron- induced inter-converting interactions and the 
produced neutron increases the ^He abundance because most of the free neutrons burned 
into "^He through D. Then PBH abundance is strongly constrained by the upper bound of 
the observational "^He abundance. For tbh ~ 10^ sec, since the produced free D can no 
longer burn into ^He, the extra free neutrons remain in D. Then (3 is severely constrained 
by the upper bound of the observational D/H. Though ''Li abundance traces D abundance 
for the longer lifetime tbh ~ 10^ sec and is produced more than SBBN at the relatively high 
7] (> 3 X 10"^''), the constraint is much weaker than D. 

In Fig. ^ we plot the most conservative bounds in the parameter region r] = (4.7 — 
8.6) X 10"^°. For tbh ~ 4 x 10^ sec such a bound is obtained by Yp at rj = 4.7 x 10~^°. For 
4 X 10^ ^ Tbh ~ 10"^ sec, the prediction of D/H at r] = 8.6 x 10~^° gives the most conservative 
bound for (3, but constraints in this region is not quantitatively accurate because we did not 
include spallation of ^He produced. 




(44) 
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Second we discuss the other case "HighD" and adopt a second set of the observational 
constraints Eqs. (ID, (H), and (H). The upper bounds for D and ^He restrict the baryon 



to photon ratio necessarily within rj = (0.90 — 8.6) x 10~^°. In Fig.^, we plot the upper 
bounds for /3 as a function of PBH mass at 77 = 2.0 x 10^^°. Because the observational 
upper bound for D is very weak, the allowed region for tbh ~ 10^ sec are enlarged. Since the 
predicted ^He abundance has the log dependence for rj, the constraint which comes from Yp 
is not changed very much compared with the case of 1] = 6.0 x 10~^°. In Fig. ^ we plot the 
most conservative upper bounds for (3 in the parameter space 1] = (0.90 — 8.6) x 10"^*^. As 
we mentioned above, the D constraint is not strict because of the weak observational upper 
bound. On the other hand, the constraint from ^Li/H is the most severe for tbh ~ 10^ sec 
because ^Li/H are almost not allowed to be produced at r] = 0.90 x 10~^°. 

Finally we estimate the number of PBHs within the neutron diffusion length dn{t) at 
t = Tbh- The neutron diffusion length is given by ^7 



dJt) = 5.8 X 10^ (^^] ' ( — ^ ' , (45) 

VlMeV/ V3 X 10-30cm; ' ^ ^ 

where at is a transport cross section. Using Eq. (|^), the number of PBHs within the horizon 
length at t = tbh is given by 

NnirBn) = (^] V = 7.7 x lO^'M^.p. (46) 

Then the number of PBHs within the neutron diffusion length is estimated by 

, 1.2 X 10^1/? , (for Tbh = 1 sec), . . 

iVd(rBH) - j 2.3 X 10^6^ , (for tbh = 10^ sec). ^^^^ 

Compared with Figs. - ^, we can see that there exist sufficiently many PBHs within the 
neutron diffusion length in the relevant parameter spaces in which we are interested and the 
number of PBHs is observationally constrained. Therefore it is justified that we can take 
the emitted hadrons as a homogeneously distributed background. 



V. CONCLUSION 

We have investigated the infiuence of hadron injection from evaporating PBHs in the 
early stage of the BBN epoch (t ~ 10~^ — 10^ sec). The above lifetimes correspond to the 
mass range of PBHs M ~ 10^ — 3 x 10^°g. When quark- ant iquark pairs or gluons are emitted 
from PBH, they immediately fragment into a lot of hadrons (pions, kaons and little nucleons) 
and the produced hadron jets are injected into the thermal plasma which is constituted by 
photons, electrons and nucleons. After the emitted hadrons are sufficiently stopped through 
the electro-magnetic interaction with the ambient photons and electrons, they can scatter 
off the ambient nucleons and they inter-convert proton and neutron each other. As a result 
more neutrons are produced and the synthesized light element abundances are drastically 
changed. In particular, ^He and D abundances are very sensitive to the neutron abundance 
and tend to increase in this scenario. Comparing with the observational data, we can 
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constrain PBH's density and their lifetime. The hadron injection which are originated by 
the direct quark- ant iquark pairs or gluons emission from PBHs has never considered at all 
in the literature. In this paper we pointed out that the hadron induced inter-converting 
interactions between protons and neutrons are very important and the energy density of 
PBH is severely constrained by the observational data. 

We obtain the following upper bounds for the initial mass fraction of PBH /3 as a function 
of the initial mass of PBH when they are formed 

P<10-'' (for 10«g < M < lO^Og) (48) 
/3<10-'' (forlO^°g<M<3xlO^°g). (49) 

Here we adopt the case of "lowD" because the produced D is more severely constrained, and 
we take the most conservative bounds which are independent of the baryon to photon ratio. 
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FIG. 1. Plot of the charged particle multiplicity {Nch) for the center of mass energy ^/s. At 
about -s/s 2± 91.2 GeV the experimental data are spread horizontally by authors. We assigned 
systematic errors of 10% for -^/s = 1.4 — 7.8 GeV . 
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FIG. 2. Upper bounds for /? which come from the observational constraints of ^He (sohd Hne), 
D (dashed hne) and ''Li (dotted hne) as a function of the PBH's mass at 7] = 6.0 X 10"^°. Here we 
take "LowD" as a deuterium observational constraints. 
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FIG. 3. The most conservative upper bounds for ''He (solid line), D (dashed line) and ''Li 
(dotted line) as mild as possible in the parameter region rjiQ = 4.7 — 8.6, where 7710 = ?] x 10^''. 
Here we take "LowD" as a deuterium observational constraint. 
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